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ABSTRACT 

The majority of cosmic rays (CRs) generated by star-forming galaxies escape those galaxies and enter the 
intergalactic medium (IGM). I show that the mean pressure of the intergalactic CRs is likely within an order 
of magnitude of the mean Lyman a forest thermal pressure. At z > 1, their pressure may have even been 
dominant. I also demonstrate that, whichever IGM phase the CRs end up in, they are a significant contributor 
to its pressure support if its temperature is ^ lO'* K, regardless of density. Which phase the CRs end up in 
depends on the structure and strength of intergalactic magnetic fields. I argue that CRs end up at least 100 kpc 
from their progenitor galaxies. CRs may self-confine in the IGM to the sound speed, generating > 10"'^ G 
magnetic fields. These considerations imply the existence and importance of a nonthermal IGM. 

Subject headings: intergalactic medium — cosmic rays — galaxies: halos 



1. INTRODUCTION 

Most of the baryonic matter in the Universe is not in stars or 
galaxies, but in the gas between the galaxies, the intergalac- 
tic medium (IGM). Far from being a uniform background 
gas, the IGM has a rich structure with several phases span- 
ning a wide range of physical conditions. About one-third of 
the baryonic mass is in the Lyman a (Lya) forest, a volume- 
filling photoionized phase with overdensities (5 < 10 and tem- 
peratures 10^ - 10^ K (Bi&Davidsen 1997). Before z~l, 
this phase made up most of the baryonic mass of the Uni- 
verse. As some of this gas collapses into smaller structures, it 
is shock heated to form the Warm Hot Intergalactic Medium 
(WHIM), a moderately dense (S - 10-100), 10^ - 10^ K 
phase that presently contains ^ 1/2 of the baryonic mass 
JCen & Ost riker 1999). Finally, some of the material has col- 
lapsed even further into "condensed" structures: these include 
the galaxies themselves, but also clusters and their hot intra- 
cluster media, as well as the multiphased gaseous haloes sur- 
rounding galaxies. The circumgalactic gas results from the 
interaction of galaxies and the IGM, through accretion onto 
galaxies and expulsion from winds, and its elf is complex, wit h 
cool dense and warm rarefied phases (e.g.. lChen et al.ll2001b . 

But is the picture of the IGM described above, an interplay 
of thermal gas, galaxies, and gravity, complete? Just as in 
galaxies, the thermal gas phases described above may be com- 
plemented by pervasive nonthermal fields: the cosmic rays 
(CRs) and magnetic fields. For the special case of galaxy clus- 
ters, detections of radio synchrotron emission conclusively 
demonstrate the presence of both cosmic rays and magnetic 
fields, and there has been much work seeking an understand- 
ing of how these phas es are generated and interact with the 
intracluster gas (e.g., ISubramanian et al.l 120061: iFerrari etld] 
|2008). Radio emission is also detected fro m galaxy filaments, 
suggesting simi lar processes ar e at work (iBrown & Rudnickl 
l2009c see also iRyu et alJ l2008i) . On a broader scale, evi- 
dence for nonthermal processes in the IGM is scant, but there 
are theoretical reasons to expect they exist. A nonthermal 
IGM can take the form of CRs accelerated by structure f or- 
mation shocks JLoeb & Waxmanll2000t[Keshet et alj|2004l) or 
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inject ed by dark mat ter annihilation throughout the Universe 
(Map elh et alJ l2006h . pervasive magnetic fields se eded by 
early galaxies or even primordially in the Big Bang (IWidrowl 
l2002h . and through the heating of the IG M by absorpt ion 
of low frequency nonthermal radio waves ( Lackill^Ol O") or 
rapid stopping of cas cades gen erated by high energy gamma 
rays (iBrodericket al.i 2012: C hang et al.ll2012l for the debate 
on whether this actually happens. see lMiniati & Elyivll2012l: 
ISchlickeiseretal .1120121) . 

A guaranteed source of nonthermal energy in the IGM is the 
CRs accelerated by star-forming galaxies. In galaxies like the 
present-day Milky Way, most of the CR energy is in protons. 
They diffuse out of the galaxy, with only a few p ercent of the 
energy lost in pion production and ionization (IStrong et al.l 
i2010i) . Starbursts and higher redshift normal galaxies have 
more gas, with stronger energy losses. However, 7-ray ob- 
servations of the starbursts M82 and NGC 253 indicate that 
^ 20-40% of their CR power is lost to pion production, 
with the majority of that power (^ 60- 80%) escaping the 
galaxy, most likely in a starburst wind (|Lacki et aP 1201 It 
lAbramowski et all l2012t lAckermann et al.l 120121) . Thus, a 
large fraction of the power injected in CRs in the Universe 
ends up in the IGM. Over the Gyrs of cosmic star formation, 
a significant intergalactic cosmic ray population grew. 

Yet the effects o f this p opulation has been largely ignored. 
iNath & BiermannI (119931) and ISamui etaH (H)05) have stud- 
ied how CRs from star-forming galaxies could have reionized 
the IGM, or heated it at high redshift. In general, though, the 
energy losses for typical CRs with kinetic energies of GeV 
or more are negligible in the present IGM. The pionic loss 
lifetime for CR protons in the z = IGM is f^r ~ 2 x 10''*(1 + 
5)~^ yr for an intergalactic density hh = 2.5(1 + (5) x 10"^ cm"^ 
jMannh eim & Schlickeiser 1994). The ionization lifetime is 
even greater for particles with kinetic energies greater than a 
few hundred MeV. A few other estimates of the intergalactic 
CR spe ctrum, with widely dispa rate assumptions, are brie fly 
given in lDar & de Rujulal (l2005f) and'Li pari & et al.l (120051) . 

CRs can nevertheless couple dynamically to the IGM and 
its weak magnetic fields. I show in this Letter that the mean 
intergalactic energy density of CRs from star-forming galax- 
ies is probably within a factor of a few of the Lya forest pres- 
sure (Section in. I also show that CRs are likely to reach the 
IGM, where they excite magnetic fluctuations in the IGM as 



they attempt to self-confine (Section[3]i. I assume Ha = 0.75, 
n„, = 0.25, //() = 70 km s"' Mpc~' for the cosmology. For 
the baryonic density, I assume fl/, = 0.045 and mean molecu- 
lar weight /i = 0.6 for a mean IGM comoving number density 
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2. THE PRESSURE OF INTERGALACTIC COSMIC RAYS 

CR nuclei experience only adiabatic losses once they leave 
a galaxy. First, CRs in a galactic wind lose momentum adi- 
abatically as the wind e xpands, sacrificin g their energy to 
push the wind outwards dVolket al.lll996l) . Galactic winds 
are common in galaxies with high specific star-formation rate, 
both high-z main sequence galaxies and compact starbursts 
(jHeckmanet al. 1990; Steidel et al. 2010). Those winds are 
powered by the efficient thermalization of supernova mechan- 
ical energy. Galaxies without supernova-launched winds may 
still have more rarefied winds originating from the halo and 
powered by CRs, as C R streaming excite plasma waves that 
push halo plasma out (iBreitschwerdt et aU Il991h . Second, 
CRs in the IGM lose momentum as the Universe expands. 

Adiabatic losses conserve the number of particles per loga- 
rithmic bin in momentum. The equation for the IGM particle 
spectrum at redshift z is therefore 
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Particles are accelerated to momenta at pinj at a rate 
dQ/d\npi^i at redshift Zinj. I relate pinj to p as p{z) = 
fadv/'injfcnjXl + z) / {I +Zinj), whcrc Eadv accounts for adiabatic 
losses within a galactic wind. I ignore ionization losses. 

CRs advected out from their host galaxy in a wind with 
density p lose momentum as p"'''-', so the Cadv factor is 
simply the cube root of the ratio of the final and ini- 
tial wind density. For example, from the Milky Way, 
the wind initially occupies the disk with radius 10 kpc 
and the CRs move with the Alfven speed 10 km s~'. If 
the wind outflows spherically and reaches 1 Mpc with a 
speed of 300 km s"' (IBreitschwerdt et al.l[l991h . then Cadv = 
[(47r(l Mpc)2(300 km s-i))/(27r(10 kpc)2(10 km s-i))]-'/^ w 
0.012. Adiabatic cooling is a much greater obstacle within 
winds from small, compact starburst regions; the original vol- 
ume is small and the termination shocks are generally far from 
the burst, resulting in a huge density contrast. 

Although winds could drain most of the CR energy by 
the time it reaches the IGM, the wind itself eventually stops 
at a termination shock when its ram pressure equals the 
surrounding IGM's pressure. The termination shock itself 
can accelerate CRs, converting eshock ^ 10% of the wind's 
kinetic energy back into CRs and releasing them into the 
IGM aokiEiL& Moi-fill 1985). For Milky Way-like galax- 
ies, the power injected at these shocks is Cshock times the lu- 
minosity in CRs in the galaxy, whereas in starburst shocks, 
the power injected is eshock times the supernova mechanical 
power, about 10 times greater per unit star formation. Reac- 
celeration can also occur at shocks within inhomogeneous 
winds ( Dorfi & Breitschwerdt2012 ). 

The injection spectrum is a power law in momentum 
(d^Q/idpinj dt) = Cp'l^), with spectral index F = 2.2. The 
normalization of the spectrum is set by a volumetric energy 
injection rate: 
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Fig. 1. — Predicted comoving cosmic CR energy density. The dasii-dotted 
line is if 10% of the original power in CRs is accelerated at termination 
shocks, while solid lines represent the original population of CRs after ex- 
periencing adiabatic losses. The grey shaded band is the expected mean IGM 
thermal pressure for T = lO" K - 10* K. 

where K is kinetic energy. I calculat e Sck by scaling to 
the c osmic star-formation rate density (iHopkins & BeacomI 
l2006h . Star formation creates young, massive stars which ac- 
celerate the CRs. The star-formation rate is directly related 
to the supernova rate as Fsn = 0.0084 yr"'(SFR/(M0 yr"')). 
Multiwavelength models of star-forming galaxies indicate 
that each supernova accelerates ^ 10^" erg of CRs (that is, 
10% of the 10^' erg released in mechanical energy). So I have 



£CR = 0.018 meV Gyr"' cm" 
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With equations [T]-[3] I calculate the evolution of the inter- 
galactic CR energy density Ucr with redshift. The solid lines 
in Figure[T]are the expected Ucr of the original CR population 
with eadv =0.001 - 1. In addition, I plot the Ucr from CRs ac- 
celerated at galactic wind termination shocks with eshock = 0.1 
of the original CR power converted back into CRs (dash- 
dotted line). 

I find that with eadv = 1, the comoving CR energy den- 
sity UQ^(z)/kB is 0.02 K cm"^ at present, and it peaked at 
z = 1.0 at 0.1 K cm"-*. For termination shocks CRs with 
Cshock = 0.1, UQ^/kg is ten times smaller For comparison, the 
expected comoving thermal IGM energy densities are 0.0075 
- 0.75 K cm"^ for T = 10'*-10'' K (grey band in Figure©. We 
see that, if eadv = 1, intergalactic CRs have an energy density 
greater than that of the Lya forest. Even for a conservative 
Cshock = 0.1, Ucr is roughly equal to the Lya forest thermal 
energy density at z « 1 and are a quarter of it at z ~ 0. Thus 
CRs may play an important role in the dynamics of the IGM. 

2.1. CR importance in other IGM phases 

CRs might not actually fill the entire volume of the Uni- 
verse, but instead may be confined to one IGM phase. These 
phases can be very overdense. But then since the same CR 
power is being squeezed into a smaller confinement volume, 
the CR energy density is likewise greater. 



What is the ratio of CR energy density to thermal energy 
density in the phase the CRs end up in? Suppose the mean 
CR energy density of the Universe is {Ucr) ■ If a fraction /) of 
the CR energy occupies an IGM phase /, and if that phase has 
a cosmic filling factor cj), then the CR energy density within 
that phase is; 

f/cR = //0(f/cR>. (4) 

The gas density in that phase is limited by mass conservation 
to: 

«;< 77-</'(niGM), (5) 

lit 

where fi, is the cosmic volume-averaged fraction of the crit- 
ical density that the IGM phase accounts for The thermal 
energy density within that phase is 11/ = (3/2)«,A:7;, giving the 
ratio 



f/cR _ nh 2{Ucr) 

r^ ■''' 
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The filling factors cancel out: if CRs end up mostly in one 
phase of the IGM, then the only relevant IGM parameters are 
its mass fraction (which is always less than 1) and its tem- 
perature. The CR energy density exceeds the thermal energy 
density as long as: 



y, ^ n,,^^ 2{Ucr) 



(7) 



At z = 0, the CRs dominate the energy density of their final 
phase if 



a- 



(8) 



where $ « 0.1 is the ratio of the actual 11^^^ and its value for 
Cadv = 1 ■ About half of the z = IGM mass has a temperature 
~ 10"* K, consisting of the Lyman a forest and condensed 
haloes, indicating a major CR contribution to the pressure. 
However, half of the IGM mass is WHIM with 10^ - 10^ K 
temperatures; if the CRs end up in WHIM, they are insignifi- 
cant. 

At higher redshift, CRs are even more important because 
psFR was much greater Furthermore, an increasingly larger 
mass fraction of the gas is in the uncollapsed 10^ K Lya for- 
est. At z = 1, I find that CRs dominate the pressure of their 
final phase if 



Ti(z= l)< 1.6 X 10' K^fi^, 



(9) 



a condition which held for most of the IGM's mass. 

Thus, CRs are likely to provide significant pressure support 
in some IGM phase at z > 1 . 

3. WHERE ARE THE COSMIC RAYS? 

3.1. Do CRs really even escape their galaxies? 

Observations of quasar absorption lines indicate that winds 
in high-z st ar-forming galaxies t ransport material out to 
> 100 kpc jHeckman et al.l 119901) . Theory likewise indi- 
cates that even Milky Wa y galaxies host CR-driven winds 
JBreitschwerdt et al]|1991l) . Ultimately, a wind with asymp- 
totic speed Voo and ram luminosity E should flow out to the 



termination shock radius Rs = \ E / {'Xtiv oaPvcyi), or 



R,p^ 1.1 Mpc 



10"*" erg sec"i 
PiGu/kB 



1/2 



1000 km s- 



0.005 K cm- 



-1/2 
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It takes roughly 1 (10) Gyr for a 1000 (100) km s"' wind to 
traverse 1 Mpc. Thus, galactic winds can transport energy 
well into the IGM. 

But suppose CRs are not advected out but simply diffuse 
from their host galaxies? CR diffusion results in a net flow 
out of the Galaxy only if the energy density outside is less than 
that inside. Since CRs do actually diffuse out of the Galaxy, 
the CR energy density in the Galactic halo is less than that 
within the disk. For a steady CR luminosity, we therefore 
have 

fMW^^_ (11) 



Vn, 



Vin 



The time that the Milky Way has been forming stars, Jmw, 
is about 10 Gyr, whereas the time that CRs stay within the 
Galactic disk is only about fi,, «30 Myr (. Connelll 119981) . 
The confinernent v olume of the Galactic disk is > 350 kpc 
(ISchlickeiseil2002h . Supposing that the CRs fill a sphere with 
radius /^out. the diffusion of CRs out of the Milky Way con- 
servatively gives us: 
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If CRs are not advected away, they must diffuse far beyond 
the star-formation and gas in the Galactic disk, well into the 
circumgalactic gas. 

3.2. How far in the IGM do they go? 

The propagation of CRs depends on the magnetic field 
strength and structure in the IGM (Adams et al. 1997). Very 
little is known about the magnetic field strength, other than 
that it less than < 1 /kG on all coherenc e lengths and < 1 nG 
on large scales dNeronov & Vovkll2010l) . The most interesting 
constraints come from gamma-ray observations. The lack of 
GeV cascade emission from blazars suggests that pair e* gen- 
erated by TeV 7-rays on their way to Earth are deflected out 
of the sightline to Earth by intergalactic magnetic fields (e.g., 
iNeronov & Vovk] j2010h . This sets a lower limit of 10"'^ G 
jDermer et a l. 2011"), but the applicability of the limits has 
been disputed (Broderick et al. 2012). 

Supposing that CRs do reach the large-scale IGM, the slow- 
est they can diffuse is Bohm diffusion - one scattering per 
gyroradius. With a mean free path of the gyroradius. 
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the distance CRs can diffuse is iBohm = \fctrl 



^Bohm = 58 Mpc 
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The maximum distance that CRs propagate into the IGM is 
set by the plasma waves that they excite while they stream. 
Within galaxies, where the magnetic energy density is com- 
parable to thermal pressure, CRs s elf-confine themselves to 
speeds less than the Alfven speed (Kul srud & Pearcd !l969). 
In the classical version of CR self-confinement, the magnetic 
fluctuations from CR streaming are much smaller than the 
mean field. However, it is thought that some nonlinear ver- 
sion of the instability can amplify magnetic fluctuations until 
they are larger than the mean field to confine CRs if necessary 
dLucek & Bell 2000). 

The sound speed is much greater than the Alfven speed 
in the IGM, though. In these conditions, the streaming 
speed of CRs is probably of order the sound speed c, = 
15 km s-i ./tJWk (Holma n et al1[T979h . This weak self- 
confinement limits the distance that intergalactic CRs travel 
to: 



Smax=Cjf =0.15 MpC 
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Note that in hot phases like the WHIM, where the CR pres- 
sure is least significant, the CRs can stream out fastest. For 
T = 10'' K, CRs can traverse about 2 Mpc, bigger than the 
typical size of a WHIM structure. Streaming CRs may escape 
quickest from those structures where they are least important. 
Even if wave generation slows down CRs in the IGM, this 
is very interesting in itself: it means that CRs excite magnetic 
fields in the IGM. The CR mean free path is at least a Larmor 
radius. If their average streaming speed over a time fstream is < 
cs, then the magnetic field must be at least B > cE/iec^t^aeam), 
where E is the particle energy and e is the particle charge. The 
characteristic magnetic field strengths are 
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lOGyr 



(16) 



on some coherence length greater than ri = 
10 pc (£'/GeV)(B/10-'^ G)"'. The existence of these 
fluctuations is co nsistent with current lim its on intergalactic 
magnetic fields (iNeronov & Vovkl 120 lOi) . The magnetic 



fluctuations that CRs can excite can have no larger than 
^ y/^nUcR on energetics grounds; for z = 0, this sets an upper 
limit of 10 nG. 



4. CONCLUSIONS 

The injection of CRs by star-forming galaxies into the IGM 
is a long-neglected kind of feedback. The energy in CRs ac- 
celerated by star-forming galaxies over the Universe's life- 
time is comparable to the thermal energy in the Lya forest. 
Although the CRs may experience strong adiabatic losses as 
they are advected away from their galaxies by winds, some of 
that energy is recovered in CRs when the winds reach their 
termination shocks. The calculated pressure of intergalactic 
CRs is therefore still within an order of magnitude of that of 
the Lya forest. If the CRs are trapped in a denser phase, their 
own density is also necessarily higher, meaning that these CRs 
are important if they end up in gas with T « 10"^ K of any 
density. CRs may be especially important at z > 1, when the 
star-formation rate density was greatest and when most of the 
IGM was relatively cool. 

If standard theories of CR streaming apply to the IGM, 
the CRs affect its magnetic structure. As they self-confine 
to roughly the sound speed, they excite magnetic fluctuations 
with B > 10"'^ G. Even so, they can penetrate through up to 
100 kpc of lO'* K gas and 1 Mpc of lO*" K gas. 

Star-forming galaxies are the not the only source of cosmic 
rays in the IGM. Active galactic nuclei and structure forma- 
tion contribute CRs too. Rather, CRs from star-forming galax- 
ies demonstrate the need to consider nonthermal processes in 
the IGM. Indeed, the importance of the nonthermal IGM is 
beginning to be recognized in other contexts. The nonthermal 
fields may prove to be just as important to the IGM as in the 
ISM. 
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